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Abstract 

Given the eventuality of neutrino and muon factories in the foreseeable future, all possible 2 — > 2 
processes involving two neutrinos, whether Dirac or Majorana ones, and two charged fermions are 
considered on the basis of the most general Lorentz invariant four-fermion effective interaction 
possible, in the limit of massless particles. Such a parametrization should enable the assessment 
of the sensitivity to physics beyond the Standard Model, including the eventual discrimination 
between the Dirac or Majorana character of neutrinos, of specific experimental beam and detector 
designs. 



1 Introduction. 



As is widely appreciated, the most general Lorentz invariant four-fermion effective parametrization of 
electroweak processes has played a central role in unravelling the basic structure and chiral properties 
of this fundamental interaction. Still to this day, the original analysis of Ref.Q is used in precision 
studies of /3-decay[0, ^, aiming at identifying at low energies tell-tale signs for physics beyond the 
Standard Model (SM). Likewise at intermediate and high energies in the purely muon and tau leptonic 
sectors, a similar parametrization ||3| has become the standard^ in terms of which to confine ever 
further parameter space, hoping to uncover a lack of overlap with that of the SM. A similar approach 
is also possible for precision studies of semi-leptonic processes which involve both the muon sector and 
the first (n, d) quark generation, for example at the intermediate energies of nuclear muon capture(5|. 

With the foreseen advent of neutrino factories and muon colliders, an analogous general analysis, 
involving in particular neutrino beams, appears to be of potential interest in the design of eventual 
experiments and detectors. For instance, the possibility of intersecting neutrino beams should not be 
dismissed, especially in the eventuality of very large intensities. Indeed, in spite of small rates, if only 
a single Uaiyb —>■ ^i^^ event for instance — as opposed to VaVh —>■ ^7^^ — were to be observed, lepton 
number violation would definitely have been established, which most likely would imply the Majorana 
character of neutrinos, one of the most pressing issues in neutrino physics today. 

This note presents such an analysis, based on the most general four-fermion effective interaction 
possible of two neutrinos and two charged fermions (whether leptons or quarks) of fixed "flavours" , or 
rather more correctly, of definite mass eigenstates, solely constrained by the requirements of Lorentz 
invariance and electric charge conservation. For instance, even though this might be realized only in 
small and peculiar classes of models beyond the SM, allowance is made for the possibility that both 
the neutrino fields and their charge conjugates couple in the effective Lagrangian density. Further- 
more, the analysis is developed separately whether for Dirac or Majorana neutrinos, with the hope to 
identify circumstances under which scattering experiments involving neutrinos could help discriminate 
between these two cases through different angular correlations for differential cross sections, given the 
high rates to be expected at neutrino factories. As is well known, the "practical Dirac-Majorana con- 
fusion theorem" Q states that within the SM, namely in the limit of massless neutrinos and {V — A) 
interactions only, these two possibilities are physically totally equivalent, and hence cannot be distin- 
guished. On the other hand, relaxing the purely {V — A) structure of the electroweak interaction by 
including at least another interaction whose chirality structure is different should suffice to evade this 
conclusion, even in the limit of massless neutrinos. 

The general classes of processes considered in this note comprise neutrino pair annihilation into 
charged leptonsQ, the inverse process of neutrino pair production through lepton annihilation, and 
finally neutrino-lepton scattering. These processes will also be considered whether either one or both 
pairs of neutrino and lepton flavoursQ, {a,b) and respectively, are identical or not. The sole 

implicit assumption is that the energy available to the reaction is both sufficiently large in order 
to justify ignoring neutrino and lepton masses, and sufficiently small in order to justify the four- 
fermion parametrization of the boson exchanges responsible for the interactions. In other words, the 
calculations are all performed in the limit of zero mass for all external neutrino and lepton mass 

""^ Henceforth, the charged fermions are referred to as leptons, even though exactly the same analysis and results apply 
to quark states, with due account then for the quark colour degree of freedom and the quark structure of the hadrons 
involved. Also, charged leptons will simply be called leptons, for short. 

^In fact, our analysis considers specific mass eigenstates for the external neutrino and lepton states, in the massless 
limit, namely when all other energy scales are much larger than the masses of these particles. By abuse of language, we 
shall refer to these mass eigenstates as "flavour" ones, following a widespread usage. 
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eigenstates. Nonetheless, effects that distinguish Majorana from Dirac neutrinos should survive in 
this massless limit. In addition to being much larger than the neutrino and lepton mass scales, the 
energy scale available to the process must also be much smaller than the energy scales associated to 
the interactions modeled by the four-fermion interactions. 

In terms of classes of processes, our analysis thus covers a wide variety of possibilities, and its 
results are presented in a manner which, it is hoped, will be found readily useful for implementation 
in numerical codes, whatever a specific model for physics beyond the SM, namely a specific set of 
effective couplings parametrizing the four-fermion interactions. In this note, no attempt is made at 
developing a systematic analysis to assess the physics potential of specific processes based on some 
particular beam and detector design, whether to look for physics beyond the SM, or discriminate 
between the Dirac or Majorana character of neutrinos. The main purpose of this work is to provide 
the general parametrization that is required for such a dedicated assessment, left for future analysis. 

The note is organized as follows. The next section provides what might be called the "kine- 
matics" of the analysis, by recalling some simple facts about Dirac and Majorana fermions. SectJ^ 
discusses the general four-fermion effective Lagrangian used in our analysis. Sects.^ to ^ then list the 
results for the three classes of processes mentioned above, with Sect.|^ only superficially illustrating 
the potential reach of these types of processes at neutrino factories. Concluding remarks are presented 
in Sect.^. 

2 A Compendium of Simple Properties 
2.1 Dirac, Weyl and Majorana spinors 

The purpose of this section is to recall a series of results relevant to Dirac, Weyl and Majorana 
quantum fermionic fields, and to specify our conventions. Since all processes are considered in the 
limit of massless neutrinos and leptons, the representation of the Clifford-Dirac algebra {7^^, 7^^} = '2^9^^ 
used througout is the chiral one, which we take to be 




(1) 



fj* [i = 1, 2, 3) being of course the usual Pauli matrices (our choice of Minkowski metric signature g^u 
is (H )). The chiral projectors {ij = it) are given by 




By definition, the charge conjugation matrix C is such that 



c-Hc = 1^ , = 75^ , c-^i'^c = -r 



C-i(7'^75)C = (7^5r 



(3) 
(4) 




with 



c 



CC^ = 1 



c^c 



(5) 



and which, in the chiral representation, is given by 




(6) 
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Given a four component Dirac spinor ip, our definition of the associated charge conjugate spinor 
is such that 

V'c = V'' = ACV^^ , (7) 

where A is some arbitrary unit phase factor, whose value would depend a priori on the choice of spinor 
field (i.e. on the neutrino or lepton flavour hereafter). This freedom in the choice of phase factor 
under charge conjugation is directly related to the "creation phase factor" of Ref.0, as shown below. 

Solutions to the free massless Dirac equation may be expanded in the helicity basis, in terms of 
the following mode representation of a Dirac quantum spinor ^£){x), 



[e~'''-''u{k,ri)b{k,r]) + e'^-''v{k,ri)d^k,ri) 



'M (27r)32|A;| 

Here, the fermionic creation and annihilation operators have the Lorentz covariant normalization 

{bik,r]),b\k',r]')] = {27rf2\k\6r,,r,'S^^Hk-k') = {dik,r]),dHk' , 
while the plane wave spinors u{k,r]) and v{k,ri) are given by, 

u{k, +) = v{k, -) 

with the Pauli bi-spinors 

X+ik) -- 







u{k, 



vik,+) 





X-Ck) 



(8) 



(9) 



(10) 



e-*^/2 cose/2 
e*^/2 sin 0/2 



X-{k) 



-e-^fl^ sin 0/2 
cos 0/2 



:ii) 



such that 



k ■ Xvik) = V Xrj{k) , Xvik)xli{k) 



1 r 



1 + r]k ■ a 



(12) 



(f and 9 being of course the usual spherical angles for the unit vector k = k/\k\ with respect to the 
axes i = 1, 2, 3, namely k = (sin cos (p, sin sin ip, cos 0). 

The value of the index r/ = it coincides with the helicity of the corresponding massless one- 
particle states, and coincides of course with the chirality of the associated quantum field. Namely, 
left- or right-handed four component Weyl spinors, with r] = — et t] = + respectively, have the 
following mode decompositions 



d^k 



-ik-x 



u{k, r])b{k, 7]) + e'''''^v{k, -ri)d\k, -rj) 



(13) 



'(oo) (27r)32|A:| 
as implied by the identification 

ijrjix) = Prji^oix) . (14) 

Hence, 6^ (A;, r/) and d^{k, rj) are the creation operators of a particle and of an antiparticule, respectively, 
both of helicity rj and momentum k. 

This identification may also be established from the chiral properties of the plane wave spinors, 
Pr,u{k,7]) = u{k,ri) , Pr,u{k,-r]) = ; Pr,v{k,r])=0 , Prjv{k,-r]) = v{k,-r]) , (15) 
u{k,r])Pr, = , u{k,-r]) Pr, =u{k,-r]) ; v{k,r]) Pr, =v{k,r]) , U(/c, -r/) = , (16) 
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as well as 

-* 1 ~h f]'y^ ~* ~* 1 — 
u{k,r])u{k,r]) = jt , v{k,ri)v{k,r]) = ft . (17) 

Likewise, their properties under charge conjugation are such that 

ClI^ {k,7]) = v{k,r]) , ClF {k,r]) = u{k,r]) ; v(k,ri) = (k,r])C , u{k,r]) = v'^ (k,r])C , (18) 

these results being specific to the helicity basis. Charge conjugates of spinors are then given by, say 
for a Dirac spinor 'tp£i{x), 

rDix)= I 7-4^ E \e-'^-''Mkv)d{k,r^)+e'^-^\v{k,^)b\k,ii)\ , (19) 
(27r)32|A;| ^ ^ 

as one should expect of course. 

Finally, let us turn to Majorana spinors. As opposed to a Dirac spinor which is comprised 
of two independent Weyl spinors of opposite chiralities, namely one of each of the two fundamental 
representations of the (covering group of the) Lorentz group, 

i)D{x)=^+{x)+^.{x) , (20) 

a Majorana spinor iPm{x) is a four component spinor which is covariant under Lorentz transformations 
but which is constructed this time from a single Weyl spinor, say of left-handed chiralityQ t] = —, and 
which is invariant under charge conjugation^ 

il^Mix) = tp-{x) + tptix) , i^lfix) = XmC^/^ = i)M{x) , (21) 

where it is now emphasized that the arbitrary phase factor Am arising in the definition of spinors 
which are self-conjugate under charge conjugation may a priori be different for each Majorana field. 
Consequently, the mode expansion of a Majorana spinor in the helicity basis is of the form. 



'I'm{x)= / ^ W"^""<k,i^)a{k,i^) + e'^-''\Mv{k,i])a\k,i^) 



(22) 



where the annihilation and creation operators a{k,r]) and a^{k,r]) obey the fermionic algebra 

{a{k,i]),a^k',ri')} = {2TTf2\k\6r,,r,'S^^\k-k') . (23) 

In terms of the quanta of the basic Weyl spinor used in the construction, we thus have the following 
correspondence^, 

aik,-) : bik,-) ■ at(fc,-) : b^{k,-) , 

a{k,+) : \Md{k,+) ; a^k,+) : Xljd^k,+) , 



^Since charge conjugation exchanges left- and right-handed chiralities, the chirality of the basic Weyl spinor used in 
this construction is irrelevant to the definition of a Majorana spinor. 

''Note that a similar definition starting from a Dirac rather than a Weyl spinor might be contemplated, then leading 
however to two independent Majorana spinors, each of which is obtained in the manner just described from a single 
distinct Weyl spinor, namely = (tpD +4'd)/\/2 and ip^^^ = —i{ipD —ip't))/\/2, in complete analogy with the real and 
imaginary parts of a single complex scalar field as well as the physical interpretation of the associated quanta as being 
particles which are or not their own antiparticles. Specifically, we have ip^j^i ~ V"- + V"- ^ > V'm'' — 4'- + V"- with 
"ip-^ = {til- -f V+)/\/2, Tp-^ = ~''P+)/V^, where ipn = 4'- Setting either ip- or ip+ to zero, the Weyl spinors 

''P-\ V'-' hence also the Majorana ones V'm\ V'm'' s-re then no longer independent, leading back to the construction 
above. 

^The complex conjugate of a complex number z is denoted z* throughout. 
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which once again shows that a^(/c, rj) is the creation operator of a particle of momentum k and helicity 
r], which furthermore is in the present case also its own antiparticle. 

Note the charge conjugation phase factor Ajv/ multiplying the creation operator contribution to 
the mode expansion of the Majorana quantum field ipMix). This phase factor corresponds exactly to 
the "creation phase factor" whose role has been emphasized already in Ref.Q on different grounds, 
and again more recently||8|. 



2.2 Differential cross sections 

All 2 — > 2 processes of interest in this note are directly considered in the center-of-mass (CM) frame 
of the reaction, with a kinematics of the form 

Pi+P2^ qi + q2 , (25) 

the quantities pi^2, Qi,2 standing of course for the four-momenta of the respective in-coming and out- 
going massless particles. Given rotation invariance, and the fact that all particles are of spin 1/2 
and of zero mass, hence of definite helicity, the sole angle of relevance is the CM scattering angle 9 
between, say, the momenta pi and qi. For all the reactions listed hereafter, the same order is used 
for the pairs {pi,P2) and {qi,q2) of the initial and final particles involved, hence leading always to the 
same interpretation for this angle 9 as being the scattering angle between the first particles in these 
two pairs of in-coming and out-going particles. 

For specific external particles of definite helicity, the differential CM cross section of all such 
processes is given by 

da 1 1 , . ,,9 da 1 1 , . ,,9 . ^ 

l-^l ' ■J7IZ^ = ^^^J^\ ■ (26) 



d^l^-^ SfGAir'^s ' dcos9 Sj 32-^3 

Here, y/s stands for the total invariant energy of the reaction, with 

s = {pi+P2? = iqi+q2f , (27) 

dilgj is the solid angle associated to the outgoing particle of normalized momentum qi = qi/\qi\, 
Sf = 2orSf = l depending on whether the two particles — including their helicity — in the final state 
are identical or not, respectively, and Ai is Feynman's scattering matrix element. Thus, it is only 
through |A4p that the differential cross section depends on|^ the scattering angle 9. Furthermore, this 
expression also shows that it is sufficient for our purposes to simply determine the amplitude Ai for 
each of the relevant processes, a single complex quantity function of 9. All our results are thus listed 
in terms of the amplitude M for each process given an arbitrary combination of helicities for the 
external states. 



3 The Effective Lagrangian 

Given a choice of external states including their helicities, Feynman's amplitude M is determined 
from the interaction Lagrangian for these particles. Assuming that the energy ^/s remains much 

^Note that this fact imphes that relative angular dependencies of cross sections are energy independent (within the 
regime to which the four-fermion parametrization applies), while of course reaction rates are directly energy dependent 
with their usual linear dependency in s. 
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smaller than any of the mass scales of the fundamental interactions at work, an effective four-fermion 
parametrization of this interaction is warranted, constrained by the sole requirements of Lorentz 
invariance and electric charge conservation. Since fermion number is not necessarily conserved in 
interactions involving neutrinos, a priori one may couple equally well the neutrino fields and their 
charge conjugates to the charged fermionic fields. For the latter, the Dirac fields that will be used 
represent the usual charged leptons (or quarks), rather than their antiparticules. It is relative to this 
choice for the charged fields that the neutrino fields and their charge conjugates are thus specified. 

With this understanding in mind, we shall consider all processes involving neutrinos (or their 
antineutrinos) of definite flavours a and b as well as leptons (or their antileptons) of flavours i and j, all 
denoted as f^, Vb, and respectively. The same notation is used for the associated spinor fields, 
except for the indication of the lepton charge. Hence, the total four-fermion effective Lagrangian that 
is considered throughout in the case of Dirac neutrino fields is of the formQ 



r . 9 



2 



8M2 
where 



Cd + C'd > (28) 



Cd = Ci+C2 + C^ + Ca, (29) 
while each of the separate contributions is given by 

Ci = S'l-^'^V^P.^Ji TjPr,,ui, + V^'^'^^'V^j^Pr^Ji Tj^^Pr,,^ + ^T^'^'^^'V^af^'' P-r,A Tja^.P^.^ , (30) 

£2 = Sf'^^'l^^Pr^Ji TjPr,,Ub + V^''^'''UlrP-vJ^ ^l^.Pv,^b + \TT''''K<y^''Pva^^ T^a^^Pr^Ub , (31) 
= S2'^''"'l7^P.,Ji I]P-vA + Vt''"'l^l''P^A I]lf.P-vA + \T^^'"'-a<y'''' P-r^Ji Tja^uP-r,A , 

(32) 

£4 = Sj-'^'l^Pr^A I]P-r,A + V^-'^'^l^J^^P.^Ji I]j^P-r,A + '"''^a'^^P,,^, i]^^.P-nA , (33) 

an implicit summation over the chiralities ija and rjh being understood of course. On the other hand, 
it is important to keep in mind that no summation over the flavour indices a and 6, nor i and j is 
implied; all four of these values are fixed at the outset, keeping open still the possibility that a and b 
might be equal or not, and likewise for i and j. 

The overall normalization factor Ag'^ /8M'^ involves a dimensionless coupling constant g as well 
as a mass scale M, while the factor 4 cancels the two 1/2 factors present in the definition of the 
chiral projection operators P±-q^ and P±ri), which appear in the effective interactions. The motivation 
for this choice of normalization is that in the specific limit of the SM, the effective interaction is 
normalized precisely in this manner with g then being the SU(2)i gauge coupling constant gi and M 
the massive gauge boson mass M\y, with in particular their tree-level relation to Fermi's constant, 
Gf/V^ = fl'|,/(8M^) (see further details below in the case of the SM). 

In the above definitions, the complex coupling coefficients {"S, F, Tj^'g's 4 parametrize the most 
general four-fermion interactions possible, including the eventuality of CP violation whenever at least 
one of these coefficients is complex. The choice for the indices rja and rjb is made such that they each 
correspond to the chirality and helicity rja or rjb of the neutrino spinor fields and associated particles 
involved in the effective coupling, with the chiralities of the leptonic fields being then determined 



''The charge exchange form of these interactions is used here, but a charge conserving one could hkewise be contem- 
plated, one being related to the other through a Fierz transformation. 
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according to the selection rules governing scalar, vector or tensor couplings, namely these chiralities are 
those of the associated neutrino for vector couplings, and opposite to it for scalar an tensor couplings. 
The situation with regards to tensor couplings is particular, in that the relation a^^^^ = ie^^'^'^ Upajl 
implies that the couplings Tj'"''"' and T^'^'^'' contribute to the effective interaction only if the chiralities 
Tja and rjh are opposite, r]a = — %, while the couplings T^"''"' and T3''"''"' contribute only if ry^ = 7]h. These 
conventions and remarks are also those that apply to the by now standard four-fermion parametrization 
used in the ^ — e sector [Q]. 

In the case of Majorana neutrino fields, a similar parametrization is of application, namely 



2 



8M2 
where 



Cm + C^ , (34) 



Compared to the definitions above, and upon using the property tpl^ = ipM characterizing Majorana 
spinors, we thus have the following correspondence between the effective coupling coefficients in the 
Majorana and Dirac cases, 

g'na,rib . gVa,r]b _j_ 5'~^<»,% _|_ gVa,-11b _j_ g~'^"-^~'^b ^gg-J 

VVa,rib . yVa,r,b _^ y-ria,Vb _^ yVa-Vb _^ y-»)a,-% ^ ^37^ 

J'ria,rib . J'Va,'nb _j_ rj^-Va^-nb _|_ J^Va-Vb _j_ rp-Va-Vb ^gg-J 

These effective Lagrangians are still not yet the most general ones possible, when either a ^ b 
or i ^ j, or both. In the above, it is implicitly assumed that the flavours a and i, on the one hand, 
and b and j on the other, couple to one another in the "current x current" representation of these 
interactions. One could still add other similar terms in which the roles of the flavours a and b, say, are 
exchanged, providing still futher interactions whenever a ^ b or i ^ j. Nonetheless, such a possibility 
may easily be included in the results hereafter, since the explicit expressions for the matrix elements 
Ai, rather than the cross sections, which are provided, are linear in the coupling coefficients. 

As a final remark, let us also note that the total neutrino fermionic number is conserved in these 
effective interactions only for couplings of type 1 and 4, {5, V, T}1'^1^'' , whereas those associated to the 
couplings of type 2 and 3, {S", V, r}2''3'''', violate that quantum number by two units. 

It is of interest to determine the effective coupling coefficients in the specific case of the elec- 
troweak Standard Model, for which the normalization factor 4g2/(8jVf^) was discussed previously 
already. Due to flavour conservation rules in that instance, different situations must be distinguished, 
depending on whether only or only Zq exchanges are involved, or both. 

Purely exchange processes arise when a = i and b = j but also a ^ b and i ^ j, in which 
case the only nonvanishing effectif coupling is 

SM : a = i ; b = j ; a^b ; i^j : V^'^ = -1 . (39) 
In the case of purely Zq neutral current processes, we have for the only nonvanishing couplings, 
SM : {a = b)^{i = j) : S^'' = sin^ 9^ , " = 1 (l - 2 sin^ Ow) , (40) 
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9w being the usual electroweak gauge mixing angle. Note that in this situation, the Lagrangians Cd 
and /Z^,, or Cm and C\,j, are equal. 

Finally, charged and neutral exchanges both contribute only when a = b = i = j , in which case 
the only nonvanishing couplings are 

SM : a = b = i = j : S^'' = sin^ % , ^f'" " + i - 2sin2 9^) . (41) 

In this case as well, the Lagrangians Cd and £^^, or Cm and C\.j, are equal. 

Any extra coupling coefficient introduced beyond these ones thus corresponds to some new 
physics beyond the Standard Model. Any particular model beyond the Standard Model predicts 
specific values for a subclass of the effective couplings parametrizing the general expression being used 
here, to which the general results to be presented hereafter may thus readily be applied. 

The remainder of the calculation proceeds straightforwardly. Given any choice of external states 
for the in-coming and out-going particles with their specific helicities, the substitution of the effective 
Lagrangian operator enables the direct evaluation of the associated matrix element M using the Fock 
algebra of the creation and annihilation operators that appear in the mode expansions of the fermion 
fields. Rather than working out the quantity |A^p through the usual trace techniques, it proves 
much more efficient to simply substitute for the explicit expressions of the u{k,r]) and v{k,r]) spinors 



solving the free Dirac equation in the helicity basis and in the chiral representation, given in Sect. 2.1 
Choosing a specific CM kinematics configuration in which only the scattering angle 9 is involved for 
the reasons of rotational invariance advocated previously, one then readily obtains a single complex 
quantity, namely simply the value for the amplitude as a function of 9. This is the procedure that 
has been applied to each of the processes, leading to the results listed herafter. 



4 Neutrino Pair Annihilation 

The first general class of processes to be considered is that of neutrino annihilations into charged 
lepton (or quark) pairs. In the Dirac case, these reactions are labelled as follows, 

{ab){ij) Dirac neutrino annihilations 



abl: 


l^a + n- 




ab2: 






ab3: 


Va + T^b- 




, ab4: 


l^a + h — 


> it + ij 


ab5: 


Va + l^b- 




, ab6: 


Va + Vb- 


> It + ^7 


ab7: 






ab8: 




>^t + ^] 



while in the Majorana case, this list reduces to 

{ab){ij) Majorana neutrino annihilations 



Mabl: Va + T^b^ i~ + q ' Mab2: Ua + Vb^ q + ■ 
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Due to identical angular momentum selection rules for all these processes, the associated matrix 
element A4 is, for all these ten processes, of the form 

M^abm) = -4s(gf^) iVi<5,, {,5,,<5-''"5-;^[^nsin2 0/2 + 2V,,.i?ii(l + cos2 0/2)] 

+6abS:>:5'T,^Cn [(1 + W) - (1 - W) cos2 6/2] 
+5-/^'5~^^Var]bDi [^12 cos2 6/2 + 25^„,^,i?i2(l + sin^ 9/2)] 
+ 6l5^^';r^ailbDiCi2 [(1 + ^?a%) - (1 - W) sin2 6/2] } 

-As(^)N2 {^abJ-^^'d-^" [^21 COs2^/2 + 25^^,^,521(1 + sin2 0/2)] 

+5ab51/!l6^r,;C2i [(1 + w) - (1 - ??a%) sin^ 6/2] 
+5-p^5-J^^riar]bD2 [A22 sin^ 0/2 + 25^^,^,522(1 + cos^ 6/2)] 

+ S'^tKj^aribD2C22 [(1 + ??a%) " (1 " W) COS^ 0/2] } , 

where, in agreement with our conventions, 9 is the scattering angle between the incoming neutrino of 
flavour a and the produced charged lepton of flavour i. The particle helicities are rja, r]i,, rji and rjj, 
respectively. Table || lists the values for the constant phase factors A'^1,2 and -Di,2 and the subsets of 
the scalar, tensor and vector effective couplings constants, in that order, which define the quantities 
^11,12,21,22, -611,12,21,22 and Cii,i2,2i,22i whether in the case of Dirac of Majorana neutrinos. 

The overall phase and sign of this amplitude is of course irrelevant physically, and is function 
of the phase convention adopted for the external |In > and |Out > states. The latter were defined by 
having the associated creation operators acting on the vacuum state |0 > in the same order as that in 
which the corresponding particles are given in the above lists of processes. For example in the case of 
the process "abl", we have thus taken 

|In>=6i(4,r/a)fo|;(fefe,%)|0> , |Out>=6j(fi,r/i)(i](fj,r/j)|0> , (43) 
in a notation that should be self-explanatory. Similarly in the case "Mab2" for instance, 

|In >= 4(4,??a)afe(4,%) |0 > , |Out >= d\{i^,7]i)b]{^j,r]j) |0 > . (44) 
Obviously, exactly all the same conventions have been used throughout this work. 

5 Neutrino Pair Production 

Although neutrino pair production processes as such pose a genuine experimental challenge for their 
detection, as opposed to processes in which they are accompanied for instance by a photon in the final 
stateQ, — > vv'^, the corresponding list of results is provided here for completeness. All 2 — > 2 
neutrino pair production processes are labelled according to the following list when both neutrinos are 
of the Dirac type, 

{ij){ab) Dirac processes 



ijl: 


4" + 


i+ - 


■^Va + Vb 


, ij2: 


4 + 


V- 


-^T^a + l'b 


ij3: 


C + 


- 


^ Va + Vb 


, ij4: 


n+ 


^7- 


Va + Vb 


ij5: 




- 


■^Va + Ub 


ij6: 




V- 


-^I'a + l'b 


ij7: 


c + 




^ i^a + h 


, ij8: 


n+ 


^7- 


l>a + h 



while in the Majorana case 
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{ij){ab) Majorana processes 



Mijl: i- + ^+ -^va + vb , Mij2: £+ + Ij i^a + i^b ■ 
For all these ten processes, the amplitude A4 is always of the following form 

A^fo)(a6) = is[^)Ni6ij {Sab5-;''^S-^''[Ansm^0/2 + 25^^,rj,Bn{l + cos^9/2)] 

-5afe5:jr-5;|^Cn [(1 + VaVb) - (1 - VaVb) cos^ 0/2] 
+5-p^5-J'^Di [Ai2 cos2 e/2 + 2V,r,,Si2(l + sin2 6/2)] 

- 5:^15^^; D,Cu [(1 + W) - (1 - VaVb) sin^ 0/2] } 
+4s (gt) iV2 {'^afc V V [^21 cos2 9/2 + 2,5,, ,,,i?2i (1 + sin2 e/2)] 

-5ab5l6l.C2i [(1 + r/a%) - (1 - w) sin2 0/2] 
+(5-'?''<5-.'"'D2 [A22 sin2 0/2 + 25,„,,,S22(1 + cos^ 0/2)] 

- 5^^'^5^^';D2C22 [(1 + w) - (1 - r?a%) cos^ 0/2] } , 

with the same conventions as previously, in particular that is the angle between the first lepton of 
flavour i and the first produced neutrino of flavour a. The different factors and coefficients appearing 
in this expression are detailed in Table whether in the case of Dirac or Majorana neutrinos. 



6 Neutrino Scattering 

Even though it would suffice in the case of neutrino scattering onto a charged lepton to give only two 
classes of processes, for instance {ai){bj) and {aj){bi), since the two other classes could be obtained 
by appropriate permutations of indices and of the coupling coefficients with their complex conjugates, 
the results for all four classes of processes are listed nonetheless, for ease of practical use, and for 
explicit check of expressions through their symmetry properties under such permutations. 

6.1 {ai){bj) neutrino scattering processes 

In the case of neutrinos of Dirac character, the list of processes is labelled according to 

{ai){bj) Dirac processes 

ail: Va + i~ + (J , ai2: Ua + if + , 

ai3: i^a + ^ i>b + (■] , ai4: Va + if ^^b + , 

ai5: i>a + Ub + £~ , ai6: Ua + if H + , 

ai7: Ua + i~ i>b + ij , ai8: Ua + if 9b + , 

while in the Majorana case 

{ai){bj) Majorana processes 
Mail: Ua + £~ Ub + £J , Mai2: i^a + if -^Vb + £j ■ 
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The general amplitude M then reads in all ten cases as follows 

^{a^m) = 4. (sf^) 6,, {6abS'^:6^^'; [An - 25,„,_,,i?n(cos2 6/2 - sin^ 9/2)] 

+5ab5-:^^5-^^Cn [1 + r?,%(cos2 6/2 - sin^ 6/2)] 
+5-/^''S-;^'^ilaVbDi [Ai2 cos2 9/2 - 2(5^,,_^,5i2(l + sin^ 9/2)] 
+ 5l6^^';wDiCi2 [(1 + W) - (1 - W) sin2 9/2)] } 
+4s (gf^) N2 [5ab5-p^5-Ji- [A21 cos2 ^/2 - 25^,,_,,i?2i(l + sin^ 0/2)] 
+5ab5l5'J^';C2i [(1 + r/a%) - (1 - w) sin^ 0/2] 
+'j;;r<5^^ry„%Z52 [^22 - 2,5^„,_^,i?22(cos2 9/2 - sin^ 6/2)] 
+ V '^'7/'' ^'^^''^2^22 [1 + r/a%(cos2 0/2 - sin^ 0/2)] } , 

being the neutrino scattering angle. The list of factors and coefficients appearing in this expression 
is detailed in Table ^ both in the Dirac and in the Majorana case. 

6.2 {aj){bi) neutrino scattering processes 

The list of processes in the Dirac case is labelled according to 



{aj){bi) Dirac processes 



i^b + it 
while in the Majorana case 



ajl: Ua + ij 

aj3: i^a + 

aj5: i>a + iJ 

aj7: Pa + C 



i^b + 

i^b + i- 
Ub + i. 



aj2: Ua + 

aj4: Ua + 

aj6: i?a + 

aj8: Ua + 



{aj){hi) Majorana processes 



Majl: va + iJ + i~ , Maj2: i/^ + q ^Vb^ C • 

The general scattering amplitude M. is of the form 

^(ai){M) = 4s (gfj^) iVi {5,fe5,7"5-/''[Ancos2 0/2-25^„,_,,i?n(l + sin2 0/2)] 

+^ab^X^'),]Cn [(1 + W) - (1 - ??a%) sin^ 0/2] 
^^X^f^wD^ [Ai2 - 25^,,_^,Si2(cos2 0/2 - sin2 0/2)] 

+ ^^:^'^^f''wB^Cr2 [1 + ??a%(cos2 0/2 - sin2 0/2)] } 
+4s (gf^) iV2 {(^afe'^;;,'5;j; [^21 - 2<5^„,_,,i?2i(cos2 0/2 - siu^ 0/2)] 
^bab^^t^'^^^-C^x [1 + r?a%(cos2 0/2 - sin^ 0/2)] 
+ V V^«^^^2 [^22 cos2 0/2 - 25^,,_^,522(1 + sin^ 0/2)] 
+ bXb^r^^wB2C22 [(1 + W) - (1 - W) sin^ 0/2)] } , 

the angle being that of the scattered neutrino. Table ^ lists the relevant factors and coefficients both 
in the Dirac and in the Majorana case. 
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6.3 {hi){aj) neutrino scattering processes 

In the Dirac case, we have the following labelling of processes 

{bi){aj) Dirac processes 



bil: 


i^b + l, - 




bi2: 






bi3: 


Ub + i, - 




bi4: 


h + if- 




bi5: 


Ub + ^i - 




bi6: 


n + if- 




bi7: 


h + ii - 




bi8: 







while in the Majorana case 



{bi){aj) Majorana processes 



Mbil: Ub + 



Mbi2: z^b + - 



The general scattering amplitude M reads 



M 



(bi)(aj) 



4s ( 



N, 5., 



1 Oij 



{5ab5-p^5-;^' [All cos2 e/2 - 25^„,_,,i?n(l + sin^ 6/2)] 
+5abS'/,:6'T^^Cii [(1 + r^aVb) - (1 - riaTib) sin^ e/2] 



(48) 



+5l5J/;;r^aVbDi [A12 - 2V,-„,5i2(cos2 9/2 - sin^ 6/2)] 
+ S-'^^S-^^riaVbDiCu [1 + r?a% (cos2 6/2 - siiP 9/2)] } 
{SabS:^':^:^; [A21 - 25^,,_^,52i(cos2 6/2 - sin2 9/2)] 
+5,65-"" 5-^'' C21 [1 + r?a% (cos2 6/2 - sin2 6/2)] 
+5-p'^8-J'>'i^aribD2 [A22 cos2 9/2 - 25^^,_^,522(1 + sin^ 9/2)] 
+ S:i:5^^^7]ar]bD2C22 [(1 + W) - (1 - r]aVb) sin2 9/2)] } , 

9 being of course the neutrino scattering angle. The factors and coefficients appearing in this repre- 
sentation are detailed in Table |5[ 

6.4 {bj){ai) neutrino scattering processes 

Processes in the Dirac case are labelled according to 

{bj){ai) Dirac processes 



bjl 
bj3 
bj5 
bj7 



Ub + q 
Vb + q 
Vb + q 
Vb+q 





, bj2: 


ffe + 


r - 
I 






, bj4: 


i>b + 


- 




^a + t. 


, bj6: 


ffe + 


- 

1 






, bjS: 


Vb + 


- 





while in the Majorana case 



{bj){ai) Majorana processes 



Mbjl: 1/6 + 1. ^ Ua + q 



Mbj2: ub + ef ^Ua+t, 
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The general scattering amplitude M is given by 

M^bj)iai) = 4. {S^) 6,, {SabS'^tS^,'; [An - 25,„,_,,i?n(cos2 6/2 - sin^ 9/2)] 

+5ab5-:^^5-^^Cn [1 + r?,%(cos2 6/2 - sin^ 6/2)] 
+6-/^''S-;^'^wDi [Ai2 cos2 9/2 - 2(5^,,_^,5i2(l + sin^ 9/2)] 
+ 5l6^^';rjamDiCi2 [(1 + W) - (1 - W) sin^ 9/2)] } 
+4s (gf^) iV2 {<)-,fe<5,-'«'(^^/'" [^21 cos2 9/2 - 25^,,_,,i?2i(l + sin^ 9/2)] 
+5ab5l5'J^';C2i [(1 + r/a%) - (1 - w) sin^ 0/2] 
+'j;;r<5^^ry„%Z52 [A22 - 2,5^„,_^,i?22(cos2 0/2 - sin^ 6/2)] 
+ V '^'7/'' ^'^^''^2^22 [1 + r/a%(cos2 0/2 - sin^ 0/2)] } , 

with being the neutrino scattering angle. Table ^ lists the relevant factors and coefficients both in 
the Dirac and Majorana cases. 



7 Exploratory Examples 

Before turning to some simple illustrative examples of the potential physics reach of these 2^2 
processes involving neutrinos, let us point out the following simple fact. For each of the six classes 
of ten processes above, when comparing the cases with Dirac or with Majorana neutrinos, one no- 
tices that Dirac neutrino processes labelled "xyn" ("x" and "y" each being one of the neutrino or 
lepton flavour symbols and n being an integer) must be in correspondence with the Majorana process 
"Mxyl" for n = 1,3, 5, 7, and with the Majorana process "Mxy2" for n = 2, 4, 6, 8. Namely, given the 
correspondence (p4|), and for a specific choice of the external particle helicities, in each case the sum 
of the corresponding four Dirac amplitudes M must coincide with the amplitude M for the Majorana 
process through the association of coupling coefficients described in (p^)-p^). It is straightforward to 
check that this correspondence is indeed obtained. 

This fact, with in particular an identical angular parametrization of the differential cross section 
for all ten processes belonging to each one of the above six 2^2 classes, also implies that a model 
independent discrimination between the Dirac of Majorana character of neutrinos would require strin- 
gent precision requirements for systematic measurements in which different helicity combinations are 
compared, in order to isolate the relevant coefficients Aap-, B^p and Cap {o(,P = 1,2), and eventually 
determine their possible relations as applicable either to the Dirac or to the Majorana case. 

However, the four-fermion parametrization used is far more general than what is usually achieved 
in any specific model beyond the SM, and as a general rule, couplings of type 2, 3 and 4, namely S2%^1, 
^2 3'2'' '^23^' arise. Under such a situation which remains quite model independent, it 

is clear that in principle there exist large classes of processes which, simply through the angular 
dependency of their differential cross sections, should enable the experimental discrimination between 
the Dirac or Majorana neutrino character. However, without any prior knowledge of the order of 
magnitude of the coupling coefficients which determine the relative strengths of different angular 
dependencies, such a discrimination in a model independent manner requires that at least one of the 
flavour pairs {ab) or (ij) be identical. Indeed, one gets different angular contributions to the Dirac or 
Majorana amplitudes M provided only at least one of the set of terms proportional either to 6ij or to 
5ab, or both, contributes to the differential cross section. Nevertheless, this still leaves open a priori 
quite some possibilities, contingent onto the properties of the neutrino beams that would become 
available in the future, in particular their flavour and helicity contents. 
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In fact, in the latter respect, any helicity content of a neutrino beam other than left-handed 
for what is thought to be a neutrino and right-handed for what is thought to an antineutrino in the 
Dirac case, depends on possible interactions beyond the SM that might contribute to the neutrino beam 
production mechanism. In any event, such a helicity "contamination" of a beam must be expected not 
to exceed, say, one percent, given present limits on neutrino helicitiesQ. Depending on the intensities 
of beams to become available at neutrino factories, this might provide an additional aspect of physical 
interest nonetheless. For the time being however, let us conservatively assume that available beams 
would only be purely left-handed for would-be Dirac neutrinos and purely right-handed for would-be 
Dirac antineutrinos. 

The experimental possibility to eventually discriminate through neutrino annihilation and scat- 
tering processes between their Dirac or Majorana character using the difference in the angular depen- 
dency of the associated cross sections, is also contingent on the strength of any new interaction beyond 
those of the SM — as is also the neutrinoless double /3-decay process for that matter[^] — , since in 
the latter model such a possibility simply evaporates in the massless limit Hence, even though the 
possibility exists in principle, its actual experimental realization hinges, on the one hand, on suffi- 
ciently intense beams at neutrino factories to allow for reasonably precise angular measurements, and 
on the other hand, on the physical existence of a new interaction different from (V — A) that couples 
sufficiently strongly to neutrinos as compared to those of the SM. Clearly, a definite assessment of the 
physics potential of such an approach to the Dirac-Majorana neutrino issue requires a systematic and 
dedicated analysis which is not attempted here, based on actual neutrino factory designs as presently 
foreseen, and the general low-energy parametrization developed here. 

Besides the potential resolution of the Dirac-Majorana neutrino issue, intense neutrino beams 
should also help turn into reality the systematic determination of the neutrino electroweak interactions, 
in a manner similar to what has been done in the leptonic (i'^//)(efe) and semi-leptonic {ud){ei'e) 
sectors |2|, Through detailed precisions measurements in different combinations of flavour and 
helicity channels, which is also part of the neutrino oscillation programmes at neutrino factories, 
it should become possible to set ever more stringent experimental bounds on the different coupling 
constants that parametrize the general effective four-fermion interaction, and search for a lack of 
overlap with those of the SM. In the same way that Refs.0, |, |] provide the expressions of all possible 
observables in terms of such coupling coefficients for the above leptonic and semi-leptonic channels, 
the results listed in this note provide those for all 2 ^ 2 processes in which two neutrinos and two 
charged fermions take part, whatever their mass eigenstates. 

In the remainder of this section, simple illustrative examples are briefly considered, whose sole 
purpose is a first exploration of some of the above issues. The emphasis here is only on the Dirac- 
Majorana neutrino issue. 

First, let us consider the elastic scattering 

+ ^ i/^ + e" , (50) 

which is thus a reaction of type "ail" belonging to the {ai){hj) class, with a = b ^ i = j. Having 
in mind for instance the left-right symmetric extensions [|ll| of the SM, based on the gauge group 
SU{2)i X SU{2)ji X U{1)b-l, let us assume for utmost simplicity that besides the Sp~ and V^'~ 
couplings of the SM, the only nonvanishing extra couplings are Si'^ and Vj^'^ . In the simplest minded 
situation where the two chiral sectors are identical in as far as is possible in all their aspects and do 

*The i/fi component is indeed dominant in neutrino beams. 
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not mix, we thus have, 



5; 



where 



sin 9w , V-\ 



-[l-2sin2 Ow] 



S 



4, 



Ml 
Ml ' 



Isiu^Ow] , (51) 



(52) 



is the physical Ught Wi to heavy W2 squared gauge boson masses ratio, with Mi ~ 80 GeV and 
M2 > 720 GeVQ, and 6w is the usual electroweak gauge mixing angle, sin^ 9w — 0.231 Q. As- 
suming then that the in-coming neutrino is necessarily left-handed, only two of the eight pos- 
sible scattering amplitudes A4 are nonvanishing due to helicity selection rules, namely those with 
{Va,r]uVb,Vj) = (~5 ~; ~) ")> (~5 +; ~5 +)• Since the angular dependency of each of these two am- 
plitudes is identical whether the neutrinos are Dirac or Majorana, the sole difference being in their 
absolute normalization, and given the difficulty in performing absolute cross section measurements, 
let us consider the summation over all processes irrespective of the particle helicities, except of course 
for that of the in-coming f^, thus corresponding to an unpolarized measurement. In the Dirac case, 
one then finds 



\M. 



(As) 



9 



8M2 



1 + 



Re 5 



4ReV, 



[1 + cos ( 



(53) 



while in the Majorana case, 



(4.)2 ( 
1 + 



9 



[4ReV- 



Re.5'---+2Rey+-+ 



4Rey-.-+2Re5+'+ 



+ 2Re5+'+]^ 

^ (1 + cos 0)2 



(54) 



Consequently, if there are indeed interactions whose chirality structure is different from those of the 
SM, processes for Dirac or Majorana neutrinos do possess different angular properties, enabling in 
principle the discrimination between the two cases through precision measurements of the angular 
dependency of the cross section, in the present case by comparing the strength of the (1 -|- cos 6)'^ term 
to its value predicted in the SM. Unfortunately, in this specific case and under the very restrictive form 
of the coefficients Si'^ and ^i^'^ considered above, limits on the possible extra interaction are already 
such that it is too weak to render any deviation observable, the relative variation in the relevant factor 
being less than a percent given the small value for 6 < (80 GeV/720 GeV)^ ~ 1.23 x 10~^. 

Scalar or tensor couplings being typically less well constrained than vector ones, let us now 
consider the possibility of an extra scalar interaction, for either of the following two elastic scattering 
reactions, 



(55) 



Both these reactions are of the "ail" type in the {ai){bj) class, with a = b^i = j in the first case, 
and a = b = i= j in the second. Within the SM, the corresponding nonvanishing couplings are thus 



5";^ ' = sin^ 9w 



v; 



±1 
4 



1 =F 2 sin^ 6w 



(56) 



where the upper (resp. lower) sign is for the {v^e) (resp. {v^^)) reaction. Assuming that S^' is the 
sole nonvanishing extra interaction, and given a left-handed in-coming v^, one finds that whether in 
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the Dirac or Majorana case the only nonvanishing ampUtudes A4 correspond to the following helicity 
combinations: {r]a,r]i;r]b, rjj) = (— , — ; — , — ), (— , — ; +, +), (— , +; — , +). Considering again the situation 
of an unpolarized measurement, the sum over these polarization states reads, 
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8M2 



4ReV, 



+ 



ReS^ 



^ {l + cos9f + ^\S^''\^{l±cose)^} , (57) 



where in the last term the upper sign corresponds to the Dirac case, and the lower sign to the 
Majorana case. Thus once again, we see that any new interaction whose chirality structure differs 
from those of the SM leads to processes in which the angular dependency discriminates between 
Dirac and Majorana neutrinos. Taking as an illustration a value {S^' \ =0.10 which is a typical 
upper-bound on such a coupling in the leptonic (e/i) sector 0|, one finds a 10% sensitivity in the 
forward-backward asymmetry. A reasonably precise measurement of the differential cross section, and 
in particular a fit to the expected distributions in either case, thus offers the prospect to resolve the 
Dirac-Majorana neutrino issue at neutrino factories. A dedicated study should hopefully confirm the 
present exploratory assessment. 

To also highlight the potentiel interest of intersecting neutrino beams, as a final example let us 
consider the following two annihilation reactions, 



(58) 



These two processes are of the type "abl" in the {ab){ij) class, with a = b^i = j in the first case, 
and a = b = i= j in the second case, and thus again with the following interactions in the SM, 



5";^ ' = sin^ 9w 



v; 



1 

±4 



1 =F 2 sin^ 9w 



(59) 



where the upper (resp. lower) sign is associated to the first (resp. second) reaction. Assuming now 
that none of the other possible coupling coefficients belongs to the classes 5*2 3' 4", V2%,'^ and T2'^^'^^^ , 
one readily finds that given left-handed initial neutrinos only, the amplitude M for these processes 
vanishes identically in the Dirac case, but not in the Majorana case, with then a specific angular 
dependency in the latter case which is function of the interactions that might contribute. Even 
though the detection of the final state products should be straightforward, the difficulty lies of course 
in the density of the initial neutrino beams even for very intense ones, implying thus an extremely low 
rate. Nonetheless, even if only through a single event, the observation of either of the above processes, 
or similar ones for other neutrino mass eigenstates, would definitely help settle the Dirac-Majorana 
neutrino puzzle through accelerator experiments. 



8 Conclusions 

In the present work, all possible 2^2 processes involving definite mass eigenstates of two neutrinos 
and two charged fermions have been considered in the massless limit, on the basis of the most general 
four-fermion effective Lagrangian possible. All interactions, whether preserving the neutrino fermion 
number or not, and for whatever helicities of the external particles, have been included. The Feynman 
amplitudes in the center-of-mass frame have been listed for all these processes, from which the relevant 
differential cross sections readily follow. Since any particular model for physics beyond the Standard 
Model predicts specific values for the four-fermion nonderivative effective coupling coefficients, this 
analysis should be of value to assess the low energy merits of any such model. 
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In the same way as has been done for /?- and //-decay and //-capture through analogous four- 
fermion effective parametrizationsQ, |5|, the advent in the foreseeable future of neutrino factories 
with their intense beams is the main motivation for the considerations developed here. Of direct 
interest is the systematic study of the electroweak interactions in the neutrino sector, by setting ever 
more stringent limits on the effective interactions of these particles through precision measurements. 
Another physics issue of great topical interest that could be addressed through such experiments is 
that of the Dirac-Majorana discrimination of the character of neutrinos. In agreement with the Dirac- 
Majorana confusion theorem|^], as soon as interactions with a chirality structure different from the 
{V — A) one of the Standard Model are introduced, there exist processes which in principle distinguish 
between these two possible characters of the neutrino through the angular dependency of differential 
cross sections, even in the massless limit. The sensitivity of such reactions however, is contingent of 
course on the relative strength of these new interactions beyond the Standard Model. Nonetheless, 
some simple examples of such a situation were briefly described, albeit not following any systematic 
investigation. 

The main purpose of this work is to provide the general results for the Feynman amplitudes for 
all possible 2^2 processes with two neutrinos. On that basis, it should now be possible to develop 
a detailed and dedicated analysis of the potential reach of different such reactions towards the above 
physics issues, given a specific design both of neutrino beams and their intensities, and of detector 
set-ups. Besides the great interest to be found in neutrino scattering experiments, the possibilities 
offered by intersecting neutrino beams should not be dismissed offhand without first a dedicated 
assessment as well, the more so since they could possibly run in parasitic mode in conjunction with 
other experiments given the proper neutrino beam geometrical layout. 
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r-i'Dh. . — fir, * 


nTih. — T)/7 
L' '/o ) '/ti 


r~i — T)h. — Tin * 

'-'2 


"^3 


C'/fa 5 '/a * 


^^Vb,Va 


B2I 


rrVb,Va* 

^3 


rpllbyVa 

-^2 


rr-Vb,Va* 


rp-Vb,Va 


rpVb,-Va* 


rpllb-Va 

^4 


n^-Vb-Va* 
^2 


^3 


J'Vb-Va* 


J'-'nb,Va 


C2I 


YVb,Va* 




ir-rib,Va* 
^4 


Y-^b,Va 


Y^b-Va* 


Y^b-'^'^ 


^2 


■fi-Vb-ria 
^3 


Y^b-Va* 


Y^'^b,Va 


D2 


1 


1 


1 






1 


1 


1 


1 


1 


A22 


gVa,Vb* 


gria,Vb 




qVa,-Vb 




<^-Va,r}b 


•^2 


•^3 




S-Va,rib 


B22 


J^Va^Vb* 


J^Va^Vb 


'-p'na,-'nt* 


rn-qa-rib 
^4 


rp-ria^Vb* 

^4 


rp-Va,rib 


^2 


-^3 


J'Va-nb* 


J'-Va,Vb 


C22 


YVa,rib* 


yVa,rib 


YVa,-Vb* 


Y^<^-^b 


Y'^i^'^b* 


Y-Va,Vb 




^3 


Y'na-Vb* 





Table 1: List of the constant factors appearing in the amphtude (^) for all {ab){ij) neutrino annihi- 
lation processes according to their labelling defined in Sect.^ 





ijl 


ij2 


ij3 


ij4 


ij5 


ij6 


ij7 


ij8 


Mijl 


Mij2 








Va 


Va 


Va 


Va 


VaK 


VaK 


VaK 


VaK 


Au 


gVb,r)a* 


gVb,Va 


g-Vb,Va* 


q-Vb,Va 


qVb-Va* 


gVb,-Va 


q-Vb,-Va* 
"^3 


q-Vb-Va 
'-'2 


g-'nb,Va* 


SVb,-Va 


Bu 


^2 


n-Vb^Va 
-^3 


rp-rib,ria* 


^4 


rpVb,-Va* 
-'4 


rpVb,-Va 


^3 


^2 


J^-ribjVa* 


J'Vb-Va 


Cn 


Y^b,Va* 


YVb,Va 


Y-Vb,ria* 


^4 


Yjlb-Va* 


Y^b-Va 


^3 


Y-'^b,~Va 


Y~Vb,Va* 


Y''^b-Va 


Di 


1 


1 




1 


1 


AaAfe 


1 


1 


1 


1 


A12 


qVa,rib* 

'-'2 


qVa,Vb 

•^3 


qVa,-Vb* 
O4 


gVa,~Vb 


g-Va,Vb* 


q-Va,Vb 

04 


•^3 


q-ria-Vb 
'-'2 


g-Va,Vb* 


gVa-rib 


B12 


rpVa,Vb* 

^2 


^3 


rpVa-Vb* 

^4 


rpria,-rib 


rp-Va,rib* 


rr-Va,Vb 

^4 


rTi-r]a,-Vb* 

-^3 


rp-Va-Vb 

^2 


J^~1la,'nb* 


'J'Va,-rib 


C12 


Y^'^'^b* 


YVa,m 


Y^a-rib* 


Y^a,-Vb 




Y-'na,Vb 


^3 








N2 






Va 


Va 


Va 


Va 


VaK 


VaK 


VaK 


VaK 


A21 


gVb,Va 


gVb,Va* 


q-Vb,Va 
04 


g-Vb,Va* 


gVb,~Va 


qVb,-Va* 


q-Vb-Va 
'-'2 


g-Vb-na* 


SVb-Va 


g-Vb,ria* 


B21 


rprit,r]a 


rprib,Va* 


^4 


rp-Vb,Va* 


rpVb-Va 


^4 


^2 


^3 


J'Vb-Va 


J'-Vb^Va* 


C21 


^3 


Y^b,Va* 


Y-^b,Va 


Y-^b,Va* 


Y^b,-Va 


jrVb-Va* 
^4 


jr-Vb-Va 
^2 


^3 


Y^b,-Va 


Y-Vb,r]a* 


D2 


1 


1 


1 


^a^l 


AflAfe 


1 


1 


1 


1 


1 


A22 


grja,Vb 


gflatVb* 


Q'na,-nb 


gVa,-Vb* 


q~Va,Vb 
O4 


g~Va,1lb* 


q-Va,-Vb 
'^2 


•^3 


gVa,-Vb 


g-Va,Vb* 


B22 


rnVa,Vb 

-^3 


^2 


rpVa,-Vb 


^4 


n^-Va,Vb 
^4 


rp-r]a,Vb* 


rp-Va,-Vb 
-'2 


rp-Va,-Vb* 

-^3 


J'Va-Vb 


J^-VatVb* 


C22 


YVa,Vb 


Y^'^'^b* 








Y-^<^'Vb* 






Y^a,~Vb 


Y-Va,rib* 



Table 2: List of the constant factors appearing in the amplitude ( ^5| ) for all {ij){ab) neutrino pair 
production processes according to their labelling defined in Sect.^. 
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ail 




dlO 


Q 1/1 


3,15 


cii 6 






iVicill 


A/Tqi9 


AT 




Va 








Vah 










An 




oVb^Va 
'-'1 


Q-VbjVa* 






•^3 


'-'1 


q-Vbj-Va 

'-'4 


/O 1 ' ILL ' 


gVb^Va- 


^11 


^4 


rp'^b,Va 
^ 1 


rp-'Hb^Va* 


rp^Vb,Va 

^2 


rpVb-Va* 

^2 


rpVb<-Va 

-^3 


1 


rp^Vb-Va 
^A 


J'-Vb-Va* 


J'VbyVa 


Cn 


^4 




•^3 


Y'Vb,Va 

2 


Y^b-Va* 

^2 


3 


^1 


•^4 


Y-^b-Va* 


Y^b,Va 




1 

X 




X 


1 


1 

X 


T 

X 




1 

X 


1 
X 


^a^b 


y4-io 


qria,Vb* 


qVa,Vb 


qVa-Vb* 

"^3 


qVa,~Vb 


q-Va,Vb* 

'~>2 


q-Va,Vb 

'-'3 


q-Va-Vb* 


q-Va,-rib 


gVa,Vb* 


g-'la-Vb 


/j 1 o 


^ 1 


^4 




rpVa-rib 


rp-Va^Vb* 

-'2 


J'-Va, rjb 


rp-'na,~'ni,* 

-^4 


^ 1 


J'ria,rib* 


J'-Va-Vb 








•^3 




•^2 


•^3 


Y-'Oa-Vb* 


Y-Va-Vb 


Y^'^'^b* 




AT" 

iV2 


Vb 


Vb 




Vb\ 




VbK 


Vb 


Vb 


Vb 


Vb 


A21 


nTlh. .Tin 


nTlh .Tin * 


r-y — 77fi .Tin 

'-'2 


"^3 


nTlh. —Tin 


ryTlh. — 77 fi * 


'/O ) '/d 


t~1 — Tlh. — Tin * 

'-'1 


gVb yVa 


g 'lb ) '/a * 


B2I 




rrVb,Va* 


rp-r)b,ria 

^2 


rr~Vb,Va* 

^3 


rrVb-Va 


rpVb,-'na* 

^2 


rp-^b-ria 
-'4 


rp~Vb-Va* 


J'VbyVa 


J'-Vb-'TIa* 


C2I 


yVb,Va 


yVb,Va* 


y-Vb,Va 


Y~rit,Va* 


Y^b,~Va 


Y^b-ria* 


^^4 




Y''^byVa 


Y-^b-Va* 


D2 




1 


1 


1 


1 


1 


1 




K^b 


1 


A22 


qVa,Vb 


gVa,Vb* 


gVa-Vb 


gVa-rjb* 


q-Va,Vb 
•^3 


q-ria,Vb* 
•^2 




q-Va-Vb* 


g-Va-Vb 




B22 


rrria,Vb 


rpVa^Vb* 


j'Va-rib 


n^Va,-Vb* 
^3 


^3 


^2 


rp-Va-Vb 


^A 


J'-Vay-Vb 


T^fja ,rib 


C22 


YVa,11b 




yVa-Vb 


YVa,-Vb* 


^3 













Table 3: List of the constant factors appearing in the amphtude (4£) for all {ai){bj) neutrino scattering 
processes according to their labelling defined in Sect. 6.1. 





ajl 


aj2 


aj3 


aj4 


aj5 


aj6 


aj7 


aj8 


Majl 


Maj2 




Vb 


Vb 


VbK 


VbK 


VbK 


VbK 


Vb 


Vb 


Vb 


Vb 


Au 


gVb,Va 


gVb,Va* 


q-Vb,Va 
•-'2 


"^3 


gVb,~Va 


gVb,-Va* 


q-Vb~Va 


g-Vb-Va* 


gVb,Va 


g-Vb-Va* 


Bu 


rpVb,Va 


^A 


^2 


rp-Vb,Va* 


^3 


^2 


^A 


rp~Vb,-ria* 




J^-'nb,-11a* 


Cn 


Y^b,Va 


Y^b,Va* 


Y-Vb,Va 


^^3 


Y^b^-Va 


Y^b,-Va* 


^A 


Y_^^b-Va* 


Y''^b'Va 


Y^'^b-Va* 


Di 


Kh 


1 


1 


1 


1 


1 


1 


KK 


KK 


1 


A12 


qVa,Vb 


gVa,Vb* 


qVa-Vb 
'->2 


•^3 


•^3 


q-ria,Vb* 
'^2 


g-Va-Vb 


q-Va,~Vb* 


g-Va,-Vb 


gVa,Vb* 


B12 


rpria,Vb 

^A 


rpVa^Vb* 


rpVa-Vb 

^2 


n^Va,—Vb* 


^3 


rn-'na,Vb* 
^2 


rp-Va-Vb 


rp-r]a,-rib* 

-'4 


p-Va-Vb 




C12 


Y^o-^^b 


Y^i^'^b* 


Y^-'-ib 


Y^'^'-'^b* 


ir-Va,11b 

^3 


^2 




^ A 


Y~'^l'^~'''lb 


Y''^<^^Vb* 


N2 


Va 


Va 


VaK 


VaK 


VaK 


VaK 


Va 


Va 


Va 


Va 


A21 


gVb,Va* 


gVb,Va 


q-Vb,Va* 
■^3 


q-Vb,Va 


gVb-Va* 


gVb-Va 


g-Vb-Va* 


q-Vb-Va 
1J4 


g-Vb-Va* 


SVb,Va 


B21 


J^VbyVa* 


rpVb ,Va 


-^3 


^2 


rprit,—r]a* 


rpVb^-Va 


rp~Vb,-Va* 


^A 


j'-Vb-'na* 


J'Vb^Va 


C21 


irVb,Va* 
^A 


Y^b,Va 


^3 


^2 


Y^b,-Va* 


^3 


Y_^''^b'-Va* 


Y^^b-Va 


Y-'^b-Va* 


Y^b,Va 


D2 


1 




1 


1 


1 


1 


KK 


1 


1 


KK 


A22 


gVa,^b* 


gVa,Vb 


gVa,-rii,* 


gVa,-Vb 


q-Va,'nb* 
'-'2 


g-Va,'nb 


q-'na,-11b* 
O4 


g-'na,-Vb 


gVa,Vb* 


g-Va-Vb 


B22 


rpVa^Vb* 


rnVa'Vb 
^A 


-^3 


n-'Va,-Vb 
^2 


rp-Va,Vb* 
-'2 


rn-Va,'rib 
-^3 


rp-ria-rib* 
-^4 


rp-Va,-Vb 


j''na,'nb* 


J'-Va-Vb 


C22 


Y^<^''^b* 


Y2°- 




Y^-'-Vb 




Y-Va^Vb 








Y-Va-Vb 



Table 4: List of the constant factors appearing in the amplitude ( |47| ) for all {aj){bi) neutrino scattering 



processes according to their labelling defined in Sect. 6. 2 
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Uli 


h,^0 
L)1Z 






DIO 


DID 


hi? 


UIO 


IVl Ull 


IVl UIZ 


AT 




Va 


VaK 


Vah 








Va 




Va 


An 


'-'1 


Qllb^Va 


Q-VbjVa* 
'-'2 


•^3 


•^3 


qVbt-Va 


'-'4 


q-Vbj-Va 

'-'1 


OTIh.TIn * 


C — Tlh. — 'tin 


^11 


^ 1 


rn'^b,Va 

-^4 


rp~Vb,Va* 

^2 


rp-'nb,Va 

-^3 


rpVb-Va* 


rpVb~Va 


rp-rib~r]a* 

^4 


rp^Vb-Va 


X'^b,Va* 


X^'^b-Va 


Cn 


^1 


•^4 


Y-rib,r)a* 
''2 


3 


yVb~Va* 


yVb,~Va 

' 2 


•^4 


^1 


yrib,Va* 


y-Vb-va 






1 


X 


1 

L 


1 
± 


1 

X 


1 






1 
X 


y4-io 


qVa,Vb* 


qVa,Vb 


qVa-Vb* 


qVa,~Vb 
'-'3 


q-Va,Vb* 

'-'3 


q-Va,Vb 

'~>2 


q-'Ha-Vb* 


q-Va,~nb 
O4 


g~Va,-Vb* 


gVa ,Vb 


/j 1 o 




rrVa,'r}b 
^ 1 


^2 


J^Va,-Vb 


rp-Va,rib* 

-'3 


'-p-Va,rib 


1 


^A 


J'-riaj-rib* 


J'Va ,Vb 








^2 


•^3 


Y-Va,rib* 
•^3 


^ 2 


T/-'?a,-»?6* 


Y-Va-Vb 

•^4 


y-Va,-Vb* 


yVa,Vb 


AT" 

iV2 


Vb 


Vb 






VbK 


Vb\ 




Vb 


Vb 


Vb 


A21 


nTih .Tin 
' to y 10, 


nTlh .Tin * 


ri — Tlh. .Tin 


'-'2 


nTlh. —Tin 


ryT)h. — Vn * 


— fJh . — Tin 

'-'1 


ry — ?1h . — 'Hn * 
'/o ) ' /a 


g-Vb-Va 


Q 'lb ) '/a * 


B2I 


rprib,ria 


rpVb,Va* 




-'2 


rrVb-Va 
^2 


^3 


rp-Vb-Va 


rp-Vb-ria* 
-'4 


J'-Vb^-Va 




C2I 


yVb,Va 


yVb,Va* 




^2 


yVb,-Va 


y^b-Va* 


y-Vb-Va 


^4 


y-Vb-Va 


yVb,Va* 


D2 


1 




1 


1 


1 


1 


K^b 


1 


1 


AaA^ 


A22 


gVa,Vb 




gVa-Vb 


gVa-nb* 


q-Va,Vb 
•^2 


•^3 


q-Va-Vb 




SVa,Vb 




B22 


rpria,Vb 


rpria,rib* 


j'Va-rib 


n^Va-Vb* 
^2 


-^2 


rn-Va^Vb* 

^3 


rp-Va-rib 
-'4 


rp-r)a,-rib* 


XVa^Vb 


J'-'Ha-rib* 


C22 


YVa,'rib 


yVa^Vb* 


yVa-Vb 


yVa-Vb* 




^3 


y^Va-Vb 


y-Va-Vb* 


yVa,Vb 


y-Va-Vb* 



Table 5: List of the constant factors appearing in the amphtude (45 
processes according to their labelhng defined in Sect.|6.3|. 



for all {bi){aj) neutrino scattering 





bjl 


bj2 


bj3 


bj4 


bj5 


bj6 


bj7 


bj8 


Mbjl 


Mbj2 




Vb 


Vb 


Vbh 


VbK 


VbK 


VbK 


Vb 


Vb 


Vb 


Vb 


All 


qVb,Va 
0^ 


gVb,Va* 


q-Vb,Va 
•^3 


"^2 


gVb,~Va 


gVb,~Va* 


g-Vb~Va 


q-Vb-Va* 


S-nb,-ria 


gVb,Va* 


Bii 


rnVb,Va 

^A 


rpVb^Va* 


^3 


rp-Vb,Va* 

^2 


^2 


-^3 


rp-Vb,~Va 


rn-Vb-Va* 

^A 


J'-Vb,~Va 


J^llb^ria* 


Cii 


yVb,Va 


yVb,Va* 


y~Vb,'na 


y-Vb,Va* 


yVb,-Va 


yrib-Va* 


y-Vb-Va 


y-Vb-Va* 


y-Vb-Va 


yVb,Va* 


Di 


1 


AaA^ 


1 


1 


1 


1 


KK 


1 


1 


KK 


A12 


gVa,Vb 


qVa,Vb* 


•^3 


qr]a,-Vb* 
'-'2 


q-Va,rib 

'-'2 


•^3 


q-Va-Vb 


g~Va~rib* 


gVa,Vb 


g~Va,~Vb* 


B12 


rpVa,Vb 


rpVa^Vb* 

^A 


-^3 


rriVa,—Vb* 
-'2 


rr-Va,'nb 

^2 


-^3 


rp-Va,-rib 
-'4 


rp-Va,-rib* 




'p-Va,-r]b* 


C12 


yVa,Vb 


yVa,Vb* 


yVa-rib 


yVa-Vb* 


y~Va,Vb 


^3 


^A 


y-Va-Vb* 


yVa,Vb 


y-Va-Vb* 


N2 


Va 


Va 


VaK 


Vah 


VaK 


VaK 


Va 


Va 


Va 


Va 


A21 


gVb,Va* 


qVb,Va 


q-Vb,Va* 
•^2 


q-Vb,Va 
•^3 


gVb-Va* 


gVb-Va 


q-Vb,-Va* 


g-Vb,-Va 


gVb,Va* 


S-Vb-Va 


B21 


J^VbjVa* 


J^Vb ,Vo. 


^2 


^3 


rpflbj—fla* 


^2 


^A 


rp-Vb-Va 


J'Vb^Va* 


J'-Vb-ria 


C21 


yVb,Va* 


^A 


y-^b,Va* 


^3 


^3 


yVb,-Va 


y-Vb-Va* 


y-Vb-Va 


yVb,Va* 


y-Vb-Va 


D2 


AaAfe 


1 


1 


1 


1 


1 


1 


KK 


KK 


1 


A22 


gVajVb* 


QVa,Vb 


gVa-rib* 


gria~rib 


•^3 


q-Va,'nb 


<^-ria,~Vb* 


q-Va,~Vb 


g-Va-Vb* 


gVa ,Vb 


B22 


^4 


rpVa,Vb 


^2 


n-'Va,-Vb 
-^3 


rn-Va,Vb* 
-'3 


^2 


rp-Va,-Vb* 


^A 


p-fla,-r]b* 


J'Va ,Vb 


C22 


y^a,Vb* 


yVa ,Vb 


yVa-rib* 


yVa,-Vb 


y-'na,Vb* 


y-Va,Vb 


y-Va,~Vb* 


y~Va,-r,b 


y-Va,-Vb* 


yVa,Vb 



Table 6: List of the constant factors appearing in the amplitude p9| ) for all {bj){ai) neutrino scattering 



processes according to their labelling defined in Sect. 6.4 
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